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We discuss a general physical mechanism for arbitrary control of the quantum states of multiple 
qubits in the symmetric Dicke subspace. The qubit-qubit coupling leads to unequal energy spacing 
in the symmetric Dicke subspace. This allows one to manipulate a prechosen transition with an 
external driving source, with other transitions remaining off-resonant. Any entangled state in the 

. . . symmetric Dicke subspace can be created from the initial ground state by tuning the driving source. 

04 , We illustrate the idea in cavity QED, but it should be applicable to other related systems. 

^ , 

O ■ PACS numbers: PACS number: 03.67.Mn, 03.65.Ud, 42.50.Dv 

(N ■ 

■ I. INTRODUCTION 

s; 

The control of quantum systems is fundamental to quantum optics as any experimental investigation of nonclassical 
. features relies on the ability to create and manipulate quantum states. Mathematically, one can design an unitary 
I— i' evolution operator to transform an initial state to the desired state. In physics, the main difficulty comes from the fact 
. that the desired evolution operator is limited by the attainable Hamiltonian. In the context of cavity QED, methods 
has been presented to force an quantized electromagnetic field localized in a cavity from an initial ground state to any 
quantum state [1-3]. In these schemes the quantum state of one or more atoms are manipulated in an controllable 
way and the coherence of the atom is transferred to the cavity field. The atoms act as the source, teaching the 
', cavity field to evolve to the desired state. 
Q^' For a system composed of multiple particles, the superpositions of product state leads to entanglement. There are 
various types of multi-particle entanglement and the characterization of entanglement has not been completed yet. 
^ ■ The control of quantum states of composite systems is a prerequisite for experimental study of entanglement properties 
^ ' [4,5]. Besides fundamental interest, the control of the time evolution of multi-particle systems is of importance for 
QQ . the implementation of quantum computers. The implementation of a quantum computational task corresponds to 
f — ■ ' the performance of an unitary transformation on the quantum register, which is composed of multiple quantum-bits 
On , (qubits) [6]. During the quantum logic operation, the qubits are generally in an entangled state. In essence, realizing 
a quantum computer is equivalent to controlling the time evolution of an N-qubit system. The Hilbert space increases 
\ exponentially as the number of qubits increases and the control of multi-qubit systems is very complex. As far as we 

■ know, no realistic mechanism has been proposed for creating an arbitrary entangled state for N-qubit systems. 

CN ■ In this paper we design a general interaction Hamiltonian which can drive an N-qubit system to evolve from an 
' ^ [ I initial product state to any superposition state in the symmetric Dicke subspace. The qubit-qubit coupling leads to 
^ unequal spacing of energy levels in the symmetric Dicke subspace. This allows one to arbitrarily manipulate a specific 
transition by using an external driving field, which is resonant with the prechosen transition but off-resonant with 
other transitions. Any symmetric entangled state can be obtained by appropriately adjusting the parameters of the 
\ external driving field. The well-known Greenberger-Horne Zeilinger (GHZ) states [7] and W states [8] are special 
classes of symmetric states. Our idea provides a possibility for the creation of general symmetric entangled states. 
The idea can be realized in realistic physical systems. 

The paper is organized as follows. In Sec. 2, we describe the interaction Hamiltonian for the manipulation of a 
prechosen transition in the symmetric Dicke subspace. In Sec. 3, we show that any symmetric entangled state can be 
created from the initial ground state by tuning the parameter of the external souce. In Sec. 4, we discuss the physical 
realization of the scheme in cavity QED and analyse decoherence effects. We also estimate the probability that the 
atomic system undergoes unwanted transitions. A summary appears in Sec. 5. 
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II. THE INTERACTION HAMILTONIAN 



We consider an N-qubit system. The qubits have two states |e) and \g) with energy- level difference wq. The 
Hamiltonian for the whole system is (assuming h — 1) 

H = Ho + Hi + H2, (1) 
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where 

Ho=uJoS,, (2) 

Hi=XS+S-, (3) 

H2 = e[e-'('^*-^)5+ + e-*^'^*-^'^-], (4) 



= lj2(\9,){9i\-\ej){ej\), (5) 

N 
N 

S~ = ^\93){ej\- 

Hi describes the qubit-qubit coupHng with couphng strength A, and H2 describes the interaction between the qubit 
system and an external driving source with frequency w. The parameters w, £, and 9 are controllable via the external 
driving source. The operators 5+, S~, and Sz obey the commutation relations of the angular-momentum operators. 
If the system is initially in a proper Dicke state \J, — J + k) [9], it would evolve within a subspace of the Dicke space 
spanned by {| J, — J) , | J, — J + 1) , | J, J)}. We consider the case that the system is initially symmetric, thus the 
Hilbert space reduces to the symmetric Dicke subspace with J = N/2, formed by the Dicke states which are symmetric 
under the permutation of any two qubits. The state \J,—J + k) is a symmetric state with k particles being in the 
state |e), i.e., 

/2j\ -1/2 

\J,-J + k)=l^j ^Pj{\ei,e2,...,ek,gk+i,9k+2,-9N)), (6) 

where {Pj} denotes the set of all distinct permutations of the qubits. In the symmetric Dicke subspace the collective 
operators 5+ and S~ act as 



S+\J,-J + k) = V(2 J - k){k + 1) I- J, -J + k + 1), (7) 
S-\J,-J + k) = ^/k{2J - k + 1)\J,-J + k - 1) 

The Hamiltonian Hi does not induce transition between different symmetric Dicke states, but shift the energy-level of 
the state | J, — J + k) by fc(2 J— fc-|-l)A. Due to the qubit-qubit coupling the energy-level spacing between \N/2, M + 1) 
and I A'^/2, M) is wq -I- 2( J — fc)A, which is depending upon the excitation number of the state | J, — J -|- k). Thus the 
spacing of energy levels in the symmetric Dicke subspace becomes unequal. The detuning between the transition 
I J, — J + k) I J, — J + k + 1) and the classical source is ojq + 2( J — fc)A — ui. 

Suppose that e <^ X and co = wq + 2{J — k)\. In this case only the transition | J, — J + k) ^ \J, —J -|- A; -|- 1) is 
resonant with the external driving source, while other transitions remains far off-resonance and can be neglected. 
Therefore, the symmetric Dicke subspace further reduces to {| J, —J -|- -|- 1) , | J, — J + k)} and the Hamiltonians Hq, 
Hi, and H2 reduce to 

Hq = 2{k - J)uJo \J, -J + k) {J, -J + k\ + (2 J -2k + 1)ljo \J, -J + k + l){J,-J + k + l\, (8) 



Hi=ak\J,-J + k){J,-J+k\+ak+i \J,-J + k + l){J,-J+k + l\, (9) 

H2 = r;fe(e-*("*-^) | J, -J + k + 1) {J, -J + k\+ e'^"'~'^ \J, -J + k) (J, -J + k + (10) 

where 



ak = k{2J - k + 1)X, 
r]k = y/{2J-k){k + l)e. 



(11) 
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In the interaction picture with respect with Hq we obtain the interaction Hamiltonian 

Hi = Hi+H2,i, (12) 

where 

H2,i = ry^e'^C^--^)^* | J, -J+k + 1) (J, -J+k\+ g-^^Cfe--^)^* | j, _ j + fc) (J, -J + k + 1|]. (13) 

Taking advantage of the unequal energy spacing induced by the qubit-qubit couphng, we can selectively manipulate 
any transition in the symmetric Dicke subspace by tuning the frequency of the external field appropritely. 

III. GENERATION OF ANY SYMMETRIC ENTANGLED STATE 

The time evolution of this system is decided by Schrodingcr's equation: 



dt 



Him))- (14) 



Perform the unitary transformation 



Then we obtain 



where 



m))=e-'"^'H\t)). (15) 
M\t)) 



dt 



= J?2,/IV'(t)), (16) 



H2,i = %(e'^ \J,-J + k + l) {J, -J + k\+ e-'^ \J,-J + k) {J, -J + k + 1|). (17) 
The interaction induces the transition 

\J,-J + k)^ cos{r]k+it)e-'°"'* \J, -J+k)- je'(^-«'=+i*) sin(%+ii) | J, -J + k + 1). (18) 
Suppose that we desire to generate the superposition state 

K 

l^d) = ^4|J,-J + fc), (19) 

where dk is a complex number, i.e., dk = \dk\ e*'^'=. Without lose of generality, we here assume that do is real, i.e., 
ifo = 0. Assume that the qubit system is initially in the state \ J,^J), i.e., all the qubits are initially in the ground 
state. We divide the time interval into K subintcirvals. The duration of the kth subinterval is ifc. During the kth 
subinterval, the frequency of the classical driving source is Wfe = + 2( J — k)X. The corresponding phase is 6k- After 
the first subinterval the system evolves to 

l^i) = cos(r/iti) I J, - J) - ie'(^i-"i*i^ sin(r?iii) \J,-J+1). (20) 
Adjust the duration ti to satisfy the following condition 

cos(?7iti) = do. (21) 

Then we obtain 

l^i) = do I J, - J) - z^l-|do|V(^^-«^*^) I J, - J + 1) . (22) 
After the second subinterval the system evolves to 



IV'2) = do I J, - J) - W 1 - |fio|'e*^^^-"^*^)[cos(r?2i2)e-'"^*= | J, - J + 1) - ie^^^^""^*^' sin(r/2i2) | J, -J + 2)] (23) 



4 



Setting 

^l-|do|'cos(??2i2) = |di|, (24) 

we have 

IV'2) = do I J, - J) - i e'[^i-«i(*i+*=)l I J, - J + 1) - ^1 - |rfo|^ - \dife'^^'+^'-"'*^-'''*'^ \J, -J + 2). (25) 
After each subinterval the highest excitation number is increased by 1. The length of the kth subinterval tk satisfy 



fc-i 

^-^Idjf cos{r]ktj) = dk- (26) 
i=o 



This leads to the final state 



K 



\^pK) = dQ\J,-J)+Y,Ck\J,-J + k), (27) 



fe=i 



where 



Ck = \dk\e'^\ (28) 

k K-k+1 fe-1 

i=i i=i j=i 

Choose the phase of the driving source appropriately so that 

fc-l K-k+l 

^Ife = Vfe - 51 + H + ^'^Z^- (29) 

Then we obtain (pk = fk and \iIjk) is just the desired state \'ipd)- In the following we propose an implementation of 
the idea with a cavity QED system, but it is not restricted in cavity QED. 

IV. PHYSICAL IMPLEMENTAION IN CAVITY QED 

Wc consider N two-level atoms interacting with a quantized cavity field and driven by a weak classical field. The 
Hamiltonian is 

H = Hf + Ha_g+Ha-c, (30) 

where 

Hf = ljqSz + ijOcO^a, (31) 

H,-q=g{a+S- +aS+), (32) 

H^-o = e[e-'^^'-^^S+ + e-^('"*-^)5-], (33) 

a+ and a arc the creation and annihilation operators for the cavity field, ujq. uJc, and ut arc the frequencies for the 
atomic transition, cavity mode, and classical field, g is coupling constant between the atoms and the cavity field, and 



£ and are the Rabi frequency and phase of the classical field. In the case 6c = ojq — ^ g\/ n +1, with n being 
the mean photon number of the cavity field, there is no energy exchange between the atomic system and the cavity. 
Then the Hamiltonian Ha-q can be replaced by the effective Hamiltonian [10] 
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N 

= Ae[53(|e,-) (e,| - \g,) {g,\)a+ a + S+ S-], (34) 

where Ac = /5c- The first and second terms describe the photon- number dependent Stark shift, and the last term 

describes the dipole coupling between the atoms induced by the cavity mode. When the cavity mode is initially in 
the vacuum state |0) it will remain in the vacuum state throughout the procedure. Then the Hamiltonian Hf reduces 
to H(j of Eq. (2) and reduces to Hi of Eq. (3). The parameters w, e, and 6 are controllable by the classical field. 

Wc now discuss the experimental implementation of the proposed scheme. In recent cavity QED experiments with 
long-living Rydbcrg atomic levels coupled to a cavity mode, the coupling constant is <? = 27r x 25kHz [11,12]. The 
atomic radiative time and photon damping time are about = 3 x 10~^s and Tc = 1 x 10~^s, respectively. Set 
N = 3, 6 = lOg, and e = g/100. Suppose that the desired state is 

m = ^(|3/2, -3/2) + |3/2, -1/2)). (35) 

Then the required atom-cavity-field interaction time ist = 7r/(4ryi) = 7r/(4-\/3£) ~ 0.29 x 10~^s. In this case the decay 
time for the superposition state \ipd) is Id = Tr/3 = 10~^s. As the cavity mode is only virtually excited the effective 
dccohcrcnce rate due to the cavity decay is K = g^ jTJP' = lOHz. The infidelity induced by the decoherence is on the 
order of t/Td + tn = 3.19 x IQ-^." 

We now consider the probability that the atomic system undergoes a transition to the state |3/2, 1/2) via the off- 
resonant coupling. The detuning between the classical field and the transition |3/2, —1/2) — )• |3/2, 1/2) is 2A = Q.2g. 
The probability that the atomic system undergoes this transition is given by 



^^-1/2^1/2 ~ smyrjl^.+XH) 0.38 x IQ-^. (36) 

With all of the above mentioned nonideal situations being considered, the total error is about 3.57 x 10~^. 



V. SUMMARY 



In conclusion, wc have discussed a physical mechanism for arbitray control of quantum states for multi-qubit in 
the symmetric Dickc subspacc. As a consequence of qubit-qubit coupling, the energy levels are not equidistant in 
the symmetric Dicke subspacc. If the driving field is tuned in resonance with a specific transition in the symmetric 
Dicke subspacc, it would be off-resonant with other transitions. This allows one to selectively manipulate a prechosen 
transition. Any symmetric entangled state can be created by adjusting the driving field. The Hamiltonian can be 
realized in physical systems with qubit-qubit coupling available. We propose an implementation of the idea in cavity 
QED. The entanglement of symmetric Dicke states is robust against the loss of particles, as demonstrated in a recent 
experiment with photons [13]. The idea opens new perspectives for research of entanglement properties of general 
symmetric multi-qubit states. The idea can also be used for preparation of entangled Dicke states for two atomic 
systems. Suppose that the first atomic system is first prepared in a superposition of Dicke states through the above 
mentioned procedure. Then this atomic system is entangled with a light field via the exchange of excitations induced 
by atom-field interaction. The two atomic systems can be prepared in an entangled Dicke state by transferring 
the excitations of the light field to the second atomic system [14]. This provides a possibility for tests of quantum 
nonlocality with two entangled atomic systems. Apart from fundamental tests of quantum theory, entangled Dicke 
states are useful for quantum communication [15]. 

This work was supported by the National Natural Science Fomidation of China under Grant No. 10674025 and the 
Doctoral Foundation of the Ministry of Education of China under Grant No. 20070386002. 
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